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ABSTRACT 

Fabrication of a complex part by automated tape placement (ATP) can require starting up a 
new tape-end in the part interior, termed “start-on-the-part.” Careful thermal management of 
the starting transient is needed to achieve uniform crystallinity and inter- laminar weld strength 
— which is the objective of this modeling effort. The transient is modeled by a Fourier-Laplace 
transform solution of the time-dependent thermal transport equation in two spatial dimensions. 
The solution is subject to a quasi- steady approximation for the speed and length of the 
consolidation head. Sample calculations are done for the Langley ATP robot applying 
PEEK/carbon fiber composite and for two upgrades in robot performance. The head starts out 
almost at rest - which meets an engineering requirement for accurate placement of the new 
tape-end. The head then rapidly accelerates until it reaches its steady state speed. This rapid 
acceleration, however, violates the quasi-steady approximation, so uniform weld strength and 
crystallinity during the starting transient are not actually achieved. The solution does give the 
elapsed time and distance from start-up to validity of the quasi-steady approximation — which 
quantifies the length of the non-uniform region. The elapsed time was always less than 0.1s 
and the elapsed distance less than 1 cm. This quantification would allow the non-uniform 
region to be either trimmed away or compensated for in the design of a part. Such 
compensation would require experiments to measure the degree of non-uniformity, because the 
solution does not provide this information. The rapid acceleration suggests that the 
consolidation roller or belt be actively synchronized to avoid abrading the tape. 
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1. INTRODUCTION 

In the in situ automated tape placement (ATP) process, carbon fiber/thermoplastic composite 
tape is spooled out and progressively welded under heat and pressure to a previous layer in 
order to build up a laminated composite panel. Often the tape is laid from edge to edge on the 
supporting tool to create a uniform panel. Then the tape ends are trimmed away to e lim inate 
the effects of the starting and stopping transients. The corresponding steady state ATP 
process has been modeled extensively, as, for example, in refs. 1-6. 

Some applications, however, require the panel to have non-uniform thickness. Such can be 
fabricated by starting a new tape-end at an interior point of the panel instead of at its edge. 

This advanced procedure is called “start-on-the-part.” It causes an initial thermal transient that 
can alter the inter-laminar weld strength and crystalline state near the new tape-end. Also, the 
in itial speed must be reduced to place the new tape-end accurately. The objective of this 
research is to maintain uniform strength and crystallinity during start-on-the-part while 
allowing for reduced initial speed. 

The modeling approach is to find a quasi-steady solution of the time-dependent thermal 
transport equation and determine the effects of the thermal transient on the effective bonding 
time, t eff , and the final crystalline state. (According to reptation theory, the inter-laminar weld 

strength after wetting occurs is proportional to t eff ' ' 4 until full strength is reached (ref. 2).) 
The speed U and length x h of the consolidation head are then varied with time to achieve a 
uniform t eff and crystalline state. 

2. TIME-DEPENDENT MODEL 

The start-on-the-part configuration being modeled is shown in Fig. 1, except that the 
consolidation head actually conforms around the new tape end to make continuous contact. 
The model uses a two-dimensional approximation that is valid when the Peclet number 
Pe » 1, where Pe is the dimensionless ratio of inertial to diffusive heat transport (ref. 3). The 
origin of the (x,y) coordinate system is located at the nip, where the time-dependent thermal 
input q(t) (in W m" 1 ) is applied. For t < 0, q = 0 , and for t > 0, q = q Q = const. The 

thermal conductivities of the supporting tool and the consolidation head are taken to be the 
same as the average thermal conductivity of the laminated composite substrate. At distances 
far away from the nip, the temperatures in the tool and in the head are maintained at 
(usually room temperature). Subject to these idealizations, the thermal transport equation is 
given by 

pc„ ^ + pC r V % = K u^ + K ^ + q„S(x)S(y)H(t) (1) 

where T(x,y,t ) is the temperature, p is the density, C p is the specific heat at constant 
pressure, K n and K 22 are the tensor thermal conductivity components, 8 (x) and <5(y) 
are Dirac delta functions, and H(t ) is the unit step function. 

The thermal input q 0 is constrained so that in the steady state (t = °°) the temperature on 
the weld interface (y = 0) decreases passively to the glass transition temperature T„ at the 
downstream end of the consolidation head x h . This constraint is given by 
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r[x ;j (°°),o,°°] = 7; 


( 2 ) 


As mentioned, U and x h are both functions of t . Dimensionless (primed) quantities can be 
defined as follows: 
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The partial derivatives then become 
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As indicated by the preceding equation, the quasi-steady approximation requires that 
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Then (1) becomes 
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This equation is solvable analytically if b is a constant; i.e., U (t')x h (t') = const . Taking the 
partial t' - derivative gives 
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where r'= dT'/dt ' . Equation (6) can now be solved by the method of Fourier and Laplace 
transforms to give 


r'(x' ,y' ,t')= ~7 e xp 
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where t* is a dummy variable of integration, 
variables by 

r(i,o,°o) = i 


The constraint (2) is given in dimensionless 


(9) 


Its application to (8) yields 

CJo exp (b)K 0 (b) ( 0) 

where K 0 is a hyperbohc Bessel function. In the quasi-steady approximation, this equation is 

assumed to hold during the transient. It provides a physical means for maintaining a constant 
b by maintaining a constant thermal input. Substituting (10) into (8) gives finally the 
constrained temperature 
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3. THERMAL TRANSIENT 

Equation (11) can be used to plot the isotherms of dimensionless temperature for various 
values of b and t ' . The nature of the thermal transient in the composite substrate during 
start-on-the-part is shown in Fig. 2 for b = 100 . (The transient in the head is a mirror image.) 
For small values of t' (Fig. 2(a)), the welding heat is confined to a small region near the nip. 
As t' increases (Fig. 2(b)), the heat penetrates more deeply into the substrate and is 
transported farther downstream. Finally, the steady state is reached (Fig. 2(c)) when the 
T'= 1 isotherm intersects the x' -axis at x'= 1 , in accordance with the constraint (9). (The 
isotherms shown in Fig. 2 are exact when U and x h are both held constant at their steady 
values during the transient. For variable U and r (l such that Ux h = const., the validity of 
the isotherms depends on the validity of the quasi- steady approximation.) 

4. QUASI-STEADY APPROXIMATION (CONTINUED) 

The condition for the validity of the quasi-steady approximation is given by (4). Besides (10), 
another steady state relationship is assumed to hold during the transient. This relationship, as 
presented in ref. 5, is 
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Here, Ax' + (T') is the positive segment of the x’ -axis that is enclosed by the 
7" -isotherm, y' max (T') is the maximum penetration of the T' -isotherm into the substrate, d 
is the thickness of the tape, and A t(T') is the total period that the dimensionless temperature 
of a material point on the weld interface equals or exceeds T' . The period Al(T') includes 
the initial heating and subsequent reheating as more plies are added. Only the positive segment 
of the x' -axis is included because only the combination of heat and pressure contributes to 
welding. 


5. UNIFORM INTER-LAMINAR WELD STRENGTH 

For ideal welding, the max im um temperature approaches but remains slightly less than the 
degradation temperature T d . Under these conditions, ref. 5 shows that t eff is directly 

proportional to At(T d ') , where T d ' is the dimensionless degradation temperature. Therefore, 
uniform inter- laminar weld strength requires that A t(T d ') = const, during the starting 
transient. A more general condition is that A t(T ') = const. ; i.e. 
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where T ' is a parametric temperature. For uniform weld strength. T ;) '= T d . Since 
b = const., U{?) oc x h {t')~ l and (13) can be solved for U it') 


or for x h (f) 


where 
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The time t is related to dimensionless time t' by (3f) 
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Finally, (11) can be used to determine y' niax (T p ' , t') and Ar'_ (T p which are then 
substituted in (14c) to give a(t') , as follows: For any t' , the value of y' as a function of 
x' is given implicitly for the T p ' -isotherm by the equation 
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The max im um penetration of the T ' -isotherm occurs at the point on the isotherm where 
dy'/dx'= 0 . Taking the partial x' -derivative of (14e) with T p held fixed and setting 
dy'/dx'=0 gives 



x' 2 +y a 




+ t*-2x’ 


dt* = 0 


(14Q 


Solving (14e) and (14f) simultaneously for a given t' gives y' max (T p ' ,t') and the 
corresponding value of x ' . (This value of x' is discarded since it is not needed.) 

For a given t ' , the value Av' + ( T p ',t ') is obtained by setting y'= 0 in (1 1) 
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and solving for the positive x' -root. 

Figure 2 shows that as t' increases, the high temperature isotherms are the first to achieve 
their steady shapes. Therefore, the starting transient is very brief for achieving uniform inter- 
laminar weld strength, which is based on the highest parametric temperature T p '=T d ' . 


6. UNIFORM CRYSTALLINE STATE 


According to ref. 6, the final crystalline state is determined by the cooling rates through the 
crystallization zone T cl to T c0 , where T c0 < T c] . Achieving a uniform crystalline state during 

start-up requires that At(T c0 ') = const. All of the formulas of section 5 then apply with 
T p '= T c0 ' . Since it takes longer for the T c0 ' -isotherm to reach its steady shape, the starting 

transient will last longer and uniform weld strength will also be achieved. A more conservative 
requirement is that A t(T ') = const. , which is implemented by setting T '= T g '= 1 . 

7. EXAMPLES 

Sample start-on-the-part calculations of U (t) and x h (t) are done for the Langley ATP robot 

applying PEEK/carbon fiber composite. Present operating conditions of the robot correspond 
to b ~ 10 . Calculations are also done for two upgrades in robot performance, which 
correspond to b = 100 and b = 1000 . These upgrades result from proportional increases in 
the steady state head speed U(°°) and consolidation length x h (°°). Increases in x^ 00 ) area 
natural consequence of making the head more conformable. 

The relevant material parameters are p = 1560 kg m 3 , C p = 1425 J kg 1 K 1 , 

K n = 6 W m 1 K 1 , K 22 = 0.72 W m 1 K 1 , T d = 823 K , T c0 = 453 K , and T g = 413 K . The 
corresponding values for the parametric temperature T p are T d '= 4.42 , T c0 '= 1.33 , and 
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T '= 1 . The fixed process parameters are = 293 K , cl = 1.27 x 10 4 m , and tape width 
w = 3.175xl0“ 2 m . 

For the first calculation (b = 10) , the input power Q 0 = wq 0 = 127 W , and the steady speed 
and length of the consolidation head are given by U (°°) = 0.0425 m s' 1 and 
x h (°°) = 1.27 xl0~ 3 m . Plots of U and x h versus t are shown in Fig. 3 for the three 
values of T ' discussed. Also plotted is the elapsed distance of the head from its starting point 

S(t) = J (/ (t*)dt * . Each of the plots is divided in two by a vertical line. The quasi-steady 
approximation (4) is valid to the right of the vertical line and invalid to the left. 

In the second calculation (b = 100) , Q 0 = 397 W , U (°°) = 0.134 m s' 1 , and 

x h (°o) = 4.02x10 ' m and the results are plotted in Fig. 4. The third calculation (b = 1000) 

has Q 0 = 1,255 W , U(°°) = 0.425 ms 1 , and x h (°°) = 0.0127 m and the results are plotted 

in Fig. 5. The elapsed times and distances from the initiation of start-up until the quasi-steady 
approximation becomes valid are summarized in Tables 1 and 2. 


Table 1. Elapsed time t (in s) from start-up to validity of the quasi-steady approx. 

b 

T p '=T d '=4.42 

T p '=T c0 '= 1.33 

if 

Oq J 

if 

10 

.0217 

.058 

.080 

100 

.0058 

.030 

.058 

1000 

.0042 

.030 

.042 

Table 2. Elapsed distance S (in m) from start-up to validity of the q 

|uasi- steady approx. 

b 

T p '=T d '=4.42 

T p '=T c0 '= 1.33 

if 

Oq J 

if 

10 

.0007 

.0018 

.0024 

100 

.0006 

.0025 

.0052 

1000 

.0013 

.0083 

.0098 


8. DISCUSSION 

At first glance, the plots seem disappointing because the underlying quasi-steady 
approximation is seen to be invalid for nearly all of the transient period. Hence, the goal of 
uniform weld strength and crystallinity during the starting transient was not achieved. The 
reason is that the transients are very steep. This steepness is advantageous, however, because 
it implies that the regions of non-uniformity are small. As shown in Table 2, the non-uniform 
regions are shorter than 1 cm for the sample calculations and sometimes much shorter. This 
table quantifies how much of the new tape end needs to be machined away to achieve 
uniformity or, alternatively, how large a non-uniform area needs to be compensated for when 
designing a laminate. Since the model does not quantify the degree of non- uniformity, 
experiments are needed that measure the variations in strength and crystallinity in the non- 
uniform region. 

Figures 3 to 5 show that the head starts moving slowly and then accelerates rapidly to its 
steady speed. This slow start is exactly what is needed for accurate placement of the new tape 
end. The rapid acceleration of the head is necessary if the region of non- uniformity is to be 
kept short. Such acceleration makes it advisable to actively synchronize the consolidation 
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roller or belt so that the tape is not abraded. These figures also show that the consolidation 
head is initially long and rapidly shrinks to its steady length. In concept, a trailing shoe could 
be lowered to exert extended pressure during the starting transient and then raised upon 
reaching the steady state. However, Table 1 shows that the duration of the starting transient is 
less that 0. 1 s for the sample calculations. This short time period appears to make the concept 
of a retractable shoe impractical. 


9. CONCLUSIONS 

The modeling objectives for start-on-the-part have been partially achieved. The time- 
dependent thermal transport equation was solved subject to a quasi- steady approximation that 
was intended to give uniform weld strength and crystallinity during the starting transient. The 
resulting solution shows that the head is initially almost at rest — in agreement with the 
engineering requirement for accurate placement of the new tape-end. The head then rapidly 
accelerates until it reaches its steady speed. 

This rapid acceleration, however, violates the quasi-steady approximation, so uniform weld 
strength and crystallinity during the starting transient are not actually achieved. The solution 
does give the elapsed time and distance from start-up to validity of the quasi-steady 
approximation — which quantifies the length of the non-uniform region. This quantification 
would allow the non-uniform region to be either trimmed away or compensated for in the 
design of a panel. Such compensation would require experiments to measure the degree of 
departure from uniformity, since the solution does not provide this information. The rapid 
acceleration suggests that the consolidation roller or belt be actively synchronized to avoid 
abrading the tape. 
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Figure 1. Start-on-the-part configuration at t = 0. The consolidation head actually conforms 
around the new tape-end to make continuous contact from x = 0 to x = x h . 





(a) t'= 0.001 

Figure 2. Isotherms of dimensionless temperature V in composite substrate for 
b = Pe/ 2 = 100 and various dimensionless times ? during start-on-the-part. 
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(b) 7y=r c0 =1.33 
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(c) T p '=T g '= 1 


Figure 4. Concluded. 
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Figure 5. Model head speed U , elapsed distance S , and consolidation length x h versus 
time t for b = Pe/ 2 = 1000 and various dimensionless parametric temperatures T '. 



Figure 5. Continued. 
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Concluded. 





